Critical and noncritical long range entanglement in Klein-Gordon fields 
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We investigate the entanglement between two spatially separated intervals in the vacuum state of 
a free ID Klein-Gordon field by means of explicit computations in the continuum limit of the linear 
harmonic chain. We demonstrate that the entanglement, which we quantify by the logarithmic 
negativity, is finite with no further need for renormalization. We find that the quantum correlations 
are scale-invariant and are determined by a function depending on the ratio of distance to length 
only. They decay much faster than the classical correlations as in the critical limit long range 
entanglement decays exponentially for separations larger than the size of the blocks, while classical 
correlations follow a power law decay. With decreasing distance of the blocks, the entanglement 
diverges as a power law in the distance. The noncritical regime manifests richer behavior, as the 
entanglement depends both on the size of the blocks and on their separation. In correspondence 
with the von Neumann entropy also long-range entanglement distinguishes critical from noncritical 
systems. 



The scaling of block entanglement in both harmonic 
and spin chains has received considerable attention re- 
cently p|. The scaling of the entanglement entropy 
S(pa) = — tr[pA log pa] of a block A has been found 
to behave in a universal way in one-dimensional criti- 
cal systems. Explicit computations for the Klein-Gordon 
massless field using density matrices have shown that the 
entanglement entropy is proportional to the area of the 
boundary 0, 0]. Using general conformal field theory 
methods it has been shown that the entropy of a block 
with size I scales as (c/3) log / for bosonic fields [HIE BE}' 
where c = 1 is the central charge in the one-dimensional 
case. This result has been verified analytically 8, 9] using 
quantum information methods both in critical bosonic 
spin chains and in critical linear harmonic chains (HC) 
[lfj, where the area-law has been proven in higher di- 
mensions in [IH [l2l ] . In noncritical chains however the 
entropy saturates for blocks larger than the correlation 
length £ ~ m _1 , where m is the energy gap, given in 
dimensionless units where h = c = 1. The von Neumann 
entropy however, requires further renormalization as it 
diverges in the continuum limit. Due to this divergence 
the mutual information / = S(pa) + S(pb) — S(pab) of 
two regions A and B has been suggested as a better mea- 
sure [131 ] since it admits a finite value in the continuum 
limit. For critical fields / can still be computed using 
conformal field theory techniques, where it scales as a 
power law with the separation between the regions, and 
decays exponentially in the massive case However, 
mutual information is not a genuine measure of entangle- 
ment as it includes both classical and quantum correla- 
tions as demonstrated by the fact that it does not vanish 
for separable states |15| . The mutual information is an 
upper bound to many other entanglement measures, as, 
for example, the distillable entanglement and the relative 
entropy of entanglement [l5| . Hence we will consider in 
the following the logarithmic negativity Et,n as a quanti- 



fier of entanglement [l6| for both pure and mixed states, 
due in part to the relative ease with which it is computed. 

We present numerical evidence that the logarithmic 
negativity admits a finite value in the critical and non 
critical field limit, making it a promising candidate for 
studying entanglement of quantum fields. We study 
the scaling of entanglement between two spatially sepa- 
rated blocks in a free one-dimensional Klein-Gordon field. 
Moreover a non-zero value of the logarithmic negativity 
implies distillability 17J in Gaussian systems. At the mo- 
ment, to the best of our knowledge, there are no methods 
from conformal field theory, which enable one to compute 
the logarithmic negativity analytically, and we thus ob- 
tain numerically its asymptotic behaviors. 

In the following we summarize our main findings before 
presenting the numerical analysis. Firstly, as there are no 
length scales in the critical field, any well-defined finite 
physical property must depend only on a single param- 
eter r = d/l, where I is the length of each of the blocks 
and d is their separation. This property is automatically 
valid both for the mutual information and the logarith- 
mic negativity both being finite in the continuum limit, 
due to the following reason. When studying the critical 
case the continuum limit is taken by taking the lattice 
constant to zero and the coupling coefficient to infinity 
in a way that the propagation velocity remains constant. 
Both the negativity and the mutual information do not 
depend on the coupling and thus the continuum limit is 
taken just by increasing the number of oscillators. Thus 
only the ratio d/l can enter as a parameter in the con- 
tinuum limit. 

Secondly, while the classical correlations between two 
sites or blocks decay as a power law in the critical regime, 
we find that quantum correlations measured by Et,n de- 
cay exponentially with the distance in this regime fol- 
lowing Eln ~ e~ !3 " r for r > 1/2, where (3 C ~ 2y2 with 
accuracy better than 1%. This result improves a previ- 
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ously found lower bound, Eln ~ e~ r for bosons [18| 
and fermions [19| . As the blocks approach each other, 
i.e. r — » 0, we find that the £xat diverges as a power law 
r~", where a = j with accuracy better than 1%. Both 
the mutual information / and the logarithmic negativity 
Eln are upper bounds to the distillable entanglement. 
Here we observe that Eln is much tighter though, since I 
scales as a power law throughout, where we obtained nu- 
merically I ~ r -0-05 jj ence classical correlations exhibit 
a power law scaling while quantum correlations exhibit 
an exponential decay. 

Thirdly, since the noncritical chain has a length scale 
proportional to the inverse mass, this system is charac- 
terized by two dimcnsionless parameters d —>■ d^ 1 and 
I — ► Z£ , such that ui — 1 , where the critical behavior 
is obtained in the limit where d — > and I — > 0. We find 
that just like the block entropy, long range entanglement 
allows us to discriminate between critical and noncritical 
fields. With respect to the blocks' size I, Elm (I) saturates 
with increasing I, while it diverges for the critical field. 
Interestingly, the saturation occurs for I > l s ~ d + 1 , in 
contrast to the block entropy for which the saturation oc- 
curs for I > l s = 1 (m — 1). With respect to the blocks' 
separation Eln{<1) ~ exp(— /3 nc (l) d 2 ) for d > I > 1. As 
d — > the entanglement in both critical and noncritical 
fields exhibits a similar behaviour since in the noncritical 
regime (finite I), Eln diverges as a power law as well. 

Entanglement between groups of discrete sites has been 
discussed before in various setups, such as the (discrete) 
Bosc-Hubbard model [2(J, spin chains [2l| and the ion 
trap [Ilj]. Our work, on the other hand, studies the be- 
havior of long-range entanglement in continuous fields. 

Let us start by describing the correspondence between 
a continuous Klein-Gordon field and the discrete chain, 
and review the computation of several quantum informa- 
tion measures. The free one-dimensional Klein-Gordon 
Hamiltonian, H — J Hdx, where 



corresponds upon discretization with a spacing a to 



H = l 
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Substituting 7Tj — > pi and (f>i q%, transforming to cir- 
cular boundary conditions and writing in dimensionless 
form we find 



1 N 

H = - E (ll + Pn - aqnqn+l) , 



(3) 



where q n and p n are canonical variables (q± = q^+i), N 
is the number of oscillators in the chain and < a < 1 is 
the coupling constant. The correlation length £ in units 
of the oscillators spacing is defined as 



1 



2(1 -a)' 



(4) 



The continuum limit of the harmonic chain corresponds 
to Eq. (TT]) in the strong coupling limit, a — ► 1, given 
that N — > oo and m = N/£ is kept constant to ensure a 
constant propagation speed. The system is critical, i.e. 
m — > 0, when N < £. 

The spectrum of Eq. ([3]) is given by 



where 9k 
express 



v k = y/\ -acos9 k , (5) 
2iTk/N and k = 0, 1, . . . , N. Then we can 
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where [afc, at] = 1, and the two-point vacuum correlation 
matrices G and H are: 



(0\qiqj\0) = g^-j^Hij = (0|^|Q) = h^-jy (7) 



Throughout this paper we consider two separated 
blocks A and B with the same size I = £L and separation 
d = £D, where L and D are the number of oscillators in 
the blocks and their separation, respectively. 

For clarity we show how the von Neumann entropy S, 
the mutual information I and the logarithmic negativity 
Eln [3 may be computed efficiently for Gaussian states 

S(A) = (/(Aj + 1/2) - f(Xj - 1/2)) , (8) 



where f(x) = xlogx and Xj are the eigenvalues of 
iGaHa with Ga being the restriction of G to a block 
A. For the logarithmic negativity we find 



£ijV = -E lo S2 (min{2Aj,l}), 



(9) 



where Xj are the eigenvalues of iGAuBHAuB, where 
Haub is obtained from Haub by time-reversal in B, i.e. 
Pb -> -Pb- 

In the continuous limit of the critical field the two- 
point correlation functions g{x\ — X2), h(x\ — X2) are 



g(x) ~ .go - log \x\ and h(x) 



1 



(10) 



The two-point correlation functions in the noncritical 
case depend on the mass m, where in the asymptotic 
limit, x >> to -1 



9{x) 



-\x\fm 



and h(x) 



-\x\fm 



(ii) 



We proceed with the presentation of the numerical re- 
sults. We examine large chains with N = 2 ■ 10 4 (and 
N = 4 ■ 10 4 oscillators in order to confirm the continuum 
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limit). We begin with the critical regime where we take 
a = 1 — 10~ 12 (deep in the critical limit). In figure [1] we 
present In Eln as a function of r. For r > 0.5 the linear 
approximation practically coincides with the computed 
values, Eln^ > 0.5) = Eq ~ e _,9cr , where the obtained 
constant is (3 C ~ 2^/2 to 1% accuracy. 

In the upper inset we observe on a log-log scale the 
power law correction to the exponential approximation. 
Assuming Eln = E\ ~ r~ a e~ /3cr , we find \n(Ei/E ) ~ 
—a lnr. For r < 0.25 (lnr < —1-4) we obtain numeri- 
cally a = | to 1% accuracy (this number was also nu- 
merically observed for critical spin systems in 23]). Note 
that a is identical to the prefactor of the entanglement 
entropy S(l) in the critical field. 

In the lower inset we confirm that Eln is scale invari- 
ant, depending only on d/l in the critical limit. We plot 
hi Eln as a function of L, the number of oscillators in 
each of the blocks, such that the ratio r = D/L is kept 
constant. The plots are given for different values of r. 
The curves are approximately constant for L sufficiently 
large to correspond to the continuum limit. We have also 
verified the same scale invariance for the mutual informa- 
tion. 
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FIG. 1: (Color online). Critical HC: ln_E_Ljv as a function of 
r = d/l. For r > 0.5 the linear approximation practically coin- 
cides with the computed values, /3 C ~ 2s/2. The dotted line is 
the the overall estimation, (Eq. 1120 . Upper inset: \a.(E\/Eo) 
as a function of lnr. Lower inset: In Eln as a function of 
L for different values of r (D, the number of oscillators that 
separate the blocks, increases with L). 

For arbitrary values of r we find 

ETn™ 1 ~ {ar~ a + f(r)) e~^\ (12) 

where f(r) ~ e~ J / r . Note that as expected f(r >1)->1 
and f(r — > 0) — > 0. (Numerically we obtain 7 ~ 3/2 
and a ~ 4/3.) The dotted line in figure [T] shows Eq. 
(|12p (on logarithmic scale), and provides a very good 
approximation. 



Let us now analyze Eq. (jTSJ) with respect to the 
blocks' size L, keeping their separation Do constant. 
First we note that the first order exponential term 
E (l) ~ exp(— f3cdo/l) has a saddle point d 2 E /dl 2 = 
at I = (3 c do/2, in which the scaling turns from exponen- 
tial at I — > to a power of 2. At I ~ (3 c do, Eq(1) already 
scales logarithmically and for I /do 3> 1, Eq(1) saturates. 
However, at this limit the power law correction becomes 
the dominant factor, where Eln ~ I 1 ' 3 - As the power 
law is obtained from the slope in a log-log plot, we show 
in figure [51 d(ln 2?£jv)/d(ln Z) as a function of InL for 
several values of the separations Dq. We also add the 
saddle points at L ~ V2Do for each of the curves, which 
indicate the power of 2. In addition, asymptotically the 
plots tend to the 1/3 power. 

12r 
10- 
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FIG. 2: (Color online). Critical HC: dQxxE LN )/d(lnL) as 
a function of L for several values of the separation Do: 
20, 100, 200, 1000, 2000. The dots indicate the saddle points 
where the power is ~ 2. Eln(L/D -> 0) ~ e L , where 
Eln (L/ D > 1) ~ L 1/3 (seen for D = 20). 
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FIG. 3: (Color online). Noncritical chain. IuEln as a func- 
tion of / for various values of do- The points Eln {I = do) fit 
a linear curve (broken red) with slope ~ —2.1. 

We now turn to investigate long-range entanglement in 
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the noncritical field, characterized by a length scale m" 1 , 
where in units of the particles' spacing in the harmonic 
chain, the length scale is £, defined in Eq. ((4|. Due to 
the existence of a length scale, entanglement has to be 
characterized by two dimensionless parameters d, I given 
in units of m -1 : d = £)/£ and I = L/£. Numerically we 
have confirmed that 

E LN (t,D,L) = E LN (£x,D/x,L/x) 

in the continuum limit, i.e. when we observe no difference 
in Eln as we simultaneously increase N and a such that 
iV/£ remains constant. The noncritical regime reduces to 
the critical one if we take d — ► and I — ► 0. 

The noncritical regime is characterized by several lim- 
its depending on the size of I and d with respect to 1 and 
with respect to each other. We expect that in correspon- 
dence with the von Neumann entropy, in the limit Z > > 1 
the scaling becomes independent of I. Interestingly, we 
observe that this is indeed true only if also I > d. In figure 
[3] we plot In Eln{1) for several values of a constant sep- 
aration do. We observe that the entanglement reaches a 
constant value £7fjy(dn), and thus distinguishes noncrit- 
ical systems from critical ones. The points Eln{1 = do), 
which are added for reference, fit a linear curve (broken 
red line). Saturation occurs at I > l s , where l s is a linear 
function of do'. l s ~ 0.75 do + !• F° r small values of do 
saturation is obtained for l s ~ 1 (m = 1 in our notation) . 
Note that each of the curves starts linearly, showing that 
the entanglement increases exponentially with I for small 
values of I. 
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FIG. 4: (Color online). Noncritical chain. IuEln as a func- 
tion of d for various values of lo. We observe exponential 
decay with a quadratic function of d (which also depends on 
lo). In the intermediate saturation regime I > d > 1 the decay 
is exponentially linear with d. Inset: power law diverge in the 
d — > limit (shown in log- log scale). The power is different 
than the 1/3 in the critical limit and in general depends on 
la. 

In the opposite limit d > I we observe exponential de- 
cay as in the critical regime, but now with different ex- 
ponent, Eln ~ exp (—f3 nc (l)d 2 \ as can be seen in fig. HJ 
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FIG. 5: (Color online). Transition from critical to noncritical 
chain. IuEln as a function of d for constant values of r = d/l 
(I increases with d). Note that all curves begin with critical 
behavior as the entanglement is approximately constant. As I 
becomes close to m _1 = 1 the noncritical behaviour emerges. 
All curves with r < 1 (d < I) coincide in correspondence with 
the saturation regime, which is independent of I. The curves 
with d > I do not coincide and characterize the regime where 
E LN ~ exp(-/3 nc (Z)d 2 ). 



where ELN{d) is shown for several values of the constant 
blocks size Iq. Note that in the d — > limit the entangle- 
ment diverges as a power law, as can be seen in the inset 
in a log-log scale. We obtain that the power is different 
from the a = i in the critical limit and in general de- 
pends on l . In addition, we observe that for l ^> d > 1 
the exponential decay becomes linear with d again, cor- 
responding to the intermediate saturation regime. 

In order to observe the transition from critical to non- 
critical behavior we plot in figure[H]lni?LAr as a function 
of d for constant ratios r = d/l. Note that all curves 
begin with a critical behavior as the entanglement is ap- 
proximately constant. As / approaches m _1 = 1 the 
noncritical behavior emerges and entanglement starts to 
decay. It can be seen that all curves with r < 1 (d < I) 
coincide at a certain point. This corresponds to the 
saturation regime, in which the entanglement decays as 
Eln ~ exp(— 2.25d), independent of I. The curves with 
d > I do not coincide and characterize the regime where 
E LN - exp (-j3 nc (l)d 2 ). 

We would like to conclude with our main results. The 
logarithmic negativity, which is a genuine measure of en- 
tanglement, is finite in the continuum limit. It is distin- 
guished from the classical correlations especially in the 
critical limit, where it decays exponentially with the sep- 
aration, while classical correlations decay as a power law. 
As the blocks approach each other, the entanglement di- 
verges as a power law, where the power seems to be equal 
to the universal prefactor of the logarithmic scaling of the 
von Neumann entropy of a large block. It would be inter- 
esting to determine analytically, whether methods from 
conformal field theory may be applied to the negativity 
and obtain a — c/3. We note that much like the en- 
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tropy of entanglement of a single block, the scaling of 
long-range entanglement allows us to discriminate criti- 
cal from noncritical behaviour. Finally, we point out that 
for the critical field both logarithmic negativity and mu- 
tual information are scale invariant and depend only on 
the ratio between the distance and length of the blocks. 

Note added. — When finalizing this paper we became 
aware of the independent work on long-range entangle- 
ment in critical spin-chains drawing similar conclusions 




Acknowledgements - This work was supported by the 
EU Integrated project (QAP) and the EPSRC QIP-IRC. 
B.R. would like to acknowledge the Israel science foun- 
dation grant no. 784/06 and German-Israeli foundation 
Grant no. 1-857. A.R. acknowledges the support of EP- 
SRC project number EP/E045049/1. M.B.P. acknowl- 
edges support by the Royal Society. 



[1] J. Eisert, M. Cramer and M.B. Plenio, larXiv:0808.3773l 
[quant-ph] . 

[2] L. Bombelli, R. Koul, J. Lee and R. Sorkin, Phys. Rev. 

D 34, 373 (1986). 
[3] M. Srednicki, Phys. Rev. Lett. 71, 666 (1993). 
[4] C. Callan and F. Wilczek, Phys. Lett. B 333, 55 (1994). 
[5] D. Rabat and M. J. Strassler, Phys. Lett. B 329, 46 

(1994). 

[6] C. Holzhey, F. Larsen and F. Wilczek, Nucl. Phys. 424, 
443 (1994). 

[7] J. S. Dowker, Class. Quant. Grav. 11, L55 (1994). 
[8] V. E. Korepin, Phys. Rev. Lett. 92, 096402 (2004). 
[9] B.-Q. Jin and V. E. Korepin, J. Stat. Physics, 116, 79 
(2003). 

[10] A. Boto and B. Reznik, Phys. Rev. A 70, 052329 (2004). 
[11] M. Cramer, J. Eisert, and M. B. Plenio, Phys. Rev. 

Lett.98, 220603 (2007). 
[12] K. Audenaert, J. Eisert, M. B. Plenio, and R. F. Werner, 

Phys. Rev. A66, 042327 (2002). 
[13] H. Casini and M. Huerta, Phys. Lett B. 600, 142 (2004). 



[14] H. Casini, Class.Quant.Grav. 24, 1293 (2007). 
[15] M.B. Plenio and S. Virmani, Quant. Inf. Comp. 7, 1 
(2007). 

[16] G. Vidal and R.F. Werner, Phys. Rev. A 65, 032314 

(2002); M. B. Plenio, Phys. Rev. Lett. 95, 090503 (2005); 

J. Eisert, PhD Thesis 2001. 
[17] M. M. Wolf G. Giedke and J. I. Cirac, Phys. Rev. 

Lett. 96, 080502 (2006). 
[18] B. Reznik, A Retzker and J. Silman, Phys. Rev. A 71, 

042104 (2005). 

[19] J. Silman and B. Reznik, Phys. Rev. A 75 052307 (2007). 
[20] U. V. Poulsen, T. Meyer and M. Lewenstein, Phys. Rev. 

A 71, 063605 (2005). 
[21] J. Keating, F. Mezzadri and M. Novaes, Phys. Rev. A 

74, 012311 (2006). 
[22] A. Retzker, J. I. Cirac, and B. Reznik, Phys. Rev. 

Lett. 94, 050504 (2005). 
[23] H. Wichterich, J. Molina- Vilaplana and S. Bose, 

larXiv:0811.T285l [quant-ph] . 



